Scaling collapse and structure functions: Identifying 
self-affinity in finite length time series. 

S. C. Chapman, B. Hnat, G. Rowlands 1 and N. W. Watkins 2 

1 Space and Astrophysics, University of Warwick, UK 
2 British Antarctic Survey (NERC), Cambridge, UK 



Manuscript submitted to 
Nonlinear Processes in Geophysics 

May 2005 



Abstract 

Empirical determination of the scaling properties and exponents of time series 
presents a formidable challenge in testing, and developing, a theoretical understand- 
ing of turbulence and other out-of-equilibrium phenomena. We discuss the special 
case of self affine time series in the context of a stochastic process. We highlight 
two complementary approaches to the differenced variable of the data: i) attempting 
a scaling collapse of the Probability Density Functions which should then be well 
described by the solution of the corresponding Fokker-Planck equation and ii) using 
structure functions to determine the scaling properties of the higher order moments. 
We consider a method of conditioning that recovers the underlying self affine scaling 
in a finite length time series, and illustrate it using a Levy flight. 
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1 Introduction 



Theories of turbulence as applied to single point measurements in a flow concern the scaling 
properties, in a statistical sense, of differenced time series, where the Taylor hypothesis is invoked 
so that the difference between measurements at some time t and a later time t + r acts as a proxy 
for the difference between measurements made at two points in the fluid separated by length scale 
L. Studies of scali ng in solar wind turbulen c e have focused on the power spectra and the structure 



functions (see e.g. 



Tu and Marsch. 



(1995); 



Horburv and Balogh (1997)) and, more recently, the 



Probability Density Function (PDF) uinat et al. 



(2002, 



2003b)). 



The statistical scaling properties of time series can in general, however, be considered in a simi- 



lar manner. There is a considerable l iterature concerning sc a ling in aurora l 



and in geomagnetic indice s (such as 



(1996); 



Voros et al 



(1998); 



Tsurutani et 



Uritskv and Pudovkin\\\99%Y 



(1990); Takalo eta 



Watkins et al. 



region 



(1993); 



(2001); 



magnetometers 



Consolini 



Kovdcs et al. 



t al. 



This is motivated in part by attempts to understand th e driven magnetospheric syste m from the per 



spective of scaling due to intrinsic processes (see e.g. 



Chapman and WatkinA J2001I) and references 



therein) and their relationship to that of the turbulent s olar wind driver. Thi s neces s itates quantita 



live comparative stud i es of s c aling in time series ( e.g. 



(2000); 



Uritskv et al. 



(2001); 



Voros et al. 



(2002); 



Takalo and Timonen ( 1998); 



Hnat et al. 



Freeman et al. 



l2003ah ). Such studies can to some 



extent fruitfully consider the low order moments, whereas a particular difficulty for comparison 
of observations with models of turbulence is t hat the i nterm ittency parameter in turbulence r(2) is 
determined by the 6 th order structure function IFriscH (11 9951) . 

More recently, studies have focussed o n the scaling properties and f unctio n al form of the PDFs 



of the differenced time series (s e e e.g. 



1 200lh : lWeieel and BakeA d2003a): 



Consolini and De Michelis ( 1998); 



Stevanova et al. 



in the case of self- similarity \Hnat et al 



(2003b); 



Sorriso-Valvo et al 



( 2003 )). This lea ds to a Fokker-Planck model 



Hnat et al. 



(2005J)). 



In this paper we describe an approach to modelling such scaling data which exploits t he data's 



self-af fine property by applying the idea of coarse graining the data iWilson 



(1979); 



Sornette 



( 20001) ). here, in the time domain. This c oarse-graining can b e achieved empiric ally, from the 



data, by a scaling collapse procedure (e.e iHnat et al. 



(2003b); 



Hnat et al. 



(2005)) (section 2), 



and, then having experimentally determined the scaling exponent, we can take the approach one 
stage further and seek to describe the data by means of a particular case of a generalise d Fokker 



Sornette 



Planck equation (GFPE, section 5). We stress here that the GFPE is here, as elsewhere (eg 

applied to a much more general class of problem than the strictly equilibrium physics for 
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which the original FP E was obtained. T he GFPE represents an alternative to the fractional Fokker- 



Planck equation (e.g. IZaslavskv I (119951) ') which is also applicable in such non-equilibrium cases. 

The critical steps in this process are then (i) establishing whether a given dataset is self affine 
and (ii) determining the scaling exponent. We highlight two important issues that arise in the 
analysis of physical datasets here. 

The first of these is that SDE models for the data, and indeed, coarse graining, deal with the 
properties of an arbitrarily large dataset. We use a well understood example of a self affine time- 
series, that of ordinary Levy motion (section 3), to show how conditioning of the data is needed 
to recover the known scaling of an arbitrarily large timeseries from one finite length. We then use 
an example of a naturally occurring timeseries, that of the AE geomagnetic index, shown previ- 
ously to exhibit self affine scaling over a range of timescales, to highlight the effectiveness, and 
the limitations, of this technique. 

The second of these is that knowledge of the scaling properties of (in principle all) the non 
zero moments is needed to capture the scaling properties of a timeseries. We again use the AE 
timeseries to illustrate this point by constructing a fractional Brownian motion ffim with the same 
second moment, but with a very different PDF. 

2 Self affine time series: concepts 

From a time series x(t) sampled at times t^, that is at evenly spaced intervals A = tk — we 
can construct a differenced time series with respect to the time increment r = sA: 

y(t,T)=x{t + T)-x(t) (1) 

so that 

x(t + T) = x(t) + y(t,r) (2) 

If we consider N successive values determined at intervals of A, that is, y{t\, A)..y(tk, A)...y(tN, A), 
their sum gives: 

x{t) = J2y(tk,A) + x (3) 
l 

where xq = x(t — NA). As N — » oo the sum ® of the y tends to the original time series x{t). 

We will make two assumptions: i) that the y(t, r) is a stochastic variable so that © can be read 
as a random walk and ii) that the y are scaling with r (to be defined next). 
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By summing adjacent pairs in the sequence, for example: 

y^(t 1 ,2A) = y(t 1 ,A) + y(t 2 ,A) (4) 

one can coarsegrain (or decimate) the time series in r. This operation gives the x(t) as a random 
walk of N/2 values of y determined at intervals of r = 2 A. We can successively coarsegrain the 
sequence an arbitrary number of times: 

x(t) = y{h,A) + y(t 2 , A) H h y(t k , A) + y(t k+1 , A) H h A) (5) 

= y^(tt, 2A) + • • • + yW(t fe , 2A) + ■ • • + y (1) (% /2 , 2A) 

= y^(h, 2 n A) + ■■■ + yW(t k , 2 n A) + ■■■ + y^\t N/2n ,2 n A) 

where this procedure is understood in the renormalization sense in that both N and n can be taken 
arbitrarily large, so that a timeseries of arbitrarily large length is considered. This procedure can 
apply to a finite sized physical system of interest provided that that system supports a large range 
of spatio- temporal scales (the smallest being A, the largest, 2 n A, n large), an example of this is 
the inertial range in fluid turbulence. 

We now consider a self affine scaling with exponent a: 

y' = 2 a y, t' = 2t, (6) 

so that 

y (n) = 2 na y ^ T = ^ A (?) 

For arbitrary r we can normalize (r = rj A) and write 

y'(t,r) = r a y(t,A) (8) 

Now if the y is a stochastic variable with self affine scaling in r, there exists a self similar PDF 
which is unchanged under the transformation ©: 

P(y'r~ a )r- a = P(y) (9) 

Importantly, the y's are not necessarily Gaussian distributed stochastic variables, but do possess 
self similarity as embodied by ©. 



\m$» for N 



This property is shared by the (a— stable) Levy flights LShlesinger et al. 
oo. The special case where the y's are both independent, identically distributed (iid) and have 
finite variance corresponds to a Brownian random walk. One can show directly from the above 
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renormalization (see for example ISornetta (12000)) that the Brownian case is just the Central Limit 



Theorem with a = 1/2 and Gaussian P{y). Here, we consider time series which possess the 
properties §tj& and ©, which may have a / 1/2 and which are time stationary solutions of a 
Fokker-Planck equation. 

An important corollary of © is of the scaling of the structure functions (and moments). The 
p th moment can be written as: 



m 



v 



/OO f-OO 
P(y)y p dy = r pa / P{y')y' p dy' (10) 
-oo J — oo 



so that 

m p ~T pa (11) 

via ©.The scaling of any of the non zero moments of a self affine time series is thus sufficient to 
determine the exponent. Importantly, all the non zero moments will share this same scaling. This 
can also be appreciated directly by writing the PDF as an expansion in the moments. If we define 
the Fourier transform of the PDF P(z) of a given time series z(t) by: 

/oo 
e ikz P(z)dz (12) 
-oo 

then it is readily shown that the p th moment is given by: 

m P = i-iT^T l*=o (13) 

where d p /dk p denotes the p th derivative with respect to k. From this it follows that the PDF can 
be expressed as an expansion in the moments: 

oo 

m r 



P {k) = j2^(ikr d4) 

Hence the PDF is defined by knowledge of all the non zero moments. 



3 Testing for self affine scaling. 

3.1 Extracting the scaling of a surrogate, a finite length Levy flight. 

We now discuss methods for testing for the property © and measuring the exponent a for a given 
finite length time series. For the purpose of illustration we consider a Levy flight of index ji = 1.8 
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which i s generated from iid random d eviates by the following algorithm for the increments (the 



y s, see 



Siesert and Friedrich. 



(2004) for details): 



fu 



sin(/ir) / cos[(l — /i)r] 



(l-M) 



(15) 



(cos(r)) ^ 

where r is a uniformly distributed random variable in the range [— ir/2, it/ 2] and v is an exponen- 
tially distributed random variable with mean 1 which is independent of r. The scaling exponent a 
from ^ and © is then related to the Levy index, /j,, by a = 1/ \i. 



ferences obtained over a wide range of r (seelMantes na and Stanley (1199 5): Hn at et al. 



fJKofthePDFofdif- 


\jinat et al. 


2003al 


b) 



for examples). This corresponds to a renormalization of the data as discussed above. We first de- 
termine the scaling exponent a from one or more of the moments via dTTT> or an estimate thereof. 
In a finite length time series, one would ideally use the scaling of the peak P(y = 0, r) (that is, 
the p = — 1 moment) with r as this is better resolved statistically than the higher order moments. 
In practice however the time series y(t, r), formed from the differences of a measured quantity, 
can as y — > be dominated by observational uncertainties. 

Figure 1 shows the scaling collapse ^ applied to a numerically generated Levy flight (fT5l 
of 10 6 increments. The curves correspond to differences at values of r = mA with m = 
[6, 10, 16, 26, 42]. Error bars denote an estimate of the expected fluctuation per bin of this his- 
togram based on Gaussian statistics (a more soph i sticate d method for estimating these for the 



Levy case may be found in 



Siegert and Friedrich. 



I 20041) 1. We see that scaling collapse can be 



verified to the precision with which the PDF is locally determined statistically. The exponent 
a = 0.544 used to achieve the scaling collapse in Figure 1 was determined empirically directly 
from an analysis of this finite length time series based on the structure functions discussed below. 

As discussed above, the scaling exponent a that successfully collapses the PDF of different r 
should emerge from th e scaling of the moment s . This is often obtained via t he generalized struc - 



ture functions (see e.g. 



Hnat et al. 



Tu and Marsch. 



(1995); 



Horburv and Balogh\&199'm: 



Hnat et al 



(2003a); 



< 2Q05|) for examples) 



S p (r) =<\y(t,r)\ p >ocr 



C(p) 



(16) 



where for self affine y(t), we have = pa (for a multifractal, ((p) is approximately quadratic 
in p). From dTH the moments will in principle share this scaling provided that the moment is non- 
zero (however in a noisy signals a moment that should vanish will be dominated by the noise). In 
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principle we can obtain a from the slopes of log- log plots of the S p versus r for any p; in practice 
this is severely limited by the finite length of the dataset. 

The ((p) for the above Levy flight obtained via (fT6b are shown plotted versus p in Figure 2. On 
such a plot we would expect a straight line C(p) ~ pa but w e see here the well known result (see 



for example 



Chechkin and Gonchar. 



l 2000h : lAfafo2o 



(2000)) that for the surrogate, the Levy time 



series of finite length, there is a turnover in scaling above p = 2 which is spurious in the sense that 
it does not reflect the self affme scaling of the infinite length timeseries. 

One way to understand this spurious bifractal scaling is that in a finite length time series the 
PDF does not have sufficient statistical resolution in the tails. Infrequently occurring large events 
in the tails will tend to dominate the higher order moments. We need to eliminate those large 
events that are poorly represented statistically without distorting the scaling properties of the time 
series. For a self affine time series an estimate of the structure functions is: 

-A 



S°= [ \y\ p P(y,r)dy^<\y\P> 
J -A 



(17) 



where the limit on the integral is proportional to the standard deviation a so that A = Qa{r), 
with some Q constant. Now cr(r) ~ V<y 2 >~ r a shares the same self affine scaling with r as 
the original timeseries y(t, r), so that if S p ~ r pa under Q then, importantly, Sp ~ T pa also. 
Provided that Q can be chosen sufficiently large to capture the dynamic range of y, and provided 
that P{y) is symmetric, (fT7b will provide a good estimate of a. This is demonstrated in figure 2 
where we also show the ((p) obtained from JT7t . 

One can thus see that once a conditioning threshold is applied, the self affine scaling of the Levy 
flight is recovered and the value of the scaling exponent is insensitive to the value of Q chosen (for 
Q sufficiently large). We obtain the value of a = 0.544 used for the scaling collapse in Figure 
1 once conditioning is applied, giving an estimate of fi = 1.84, consistent with the index used 
to generate the synthetic Levy flight (fT3l . Similar results for a surrogate Levy dataset have been 
obtained by M. Parkinson (private communication, 2004). 

An analogous procedure to dl7t can also be re a lized by means of a trun cated wavelet expansion 



of the data (see for example 



Kovdcs et al. 



(200 U): \Mansenev et al. 



(200J)). 



In (fT7b we assumed self affine scaling in choosing the functional form of the limits of the 
integral. In a given time series the scaling may not be known a priori. If for example the time 
series were multifractal (((p) quadratic in p) we would obtain from (fT7b a ((p) which varied 
systematically with Q. In practice, several other factors may also be present in a time series which 
may additionally reduce the accuracy of the approximation ST7\ . 
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3.2 Extracting the scaling of a 'natural' example, the AE timeseries 



To illustrate the abo ve, we consider an i nterval of the AE index shown previously to exhibit weakly 



multifractal scaling iHnat et al. 



(2005)). The scaling inde x is not with i n the 



it has been modelled with a GFPE rather than a Levy walk 



The PDF of differenced AE is asymmetric 



Hnat 



et al. 



Levy range and thus 



(2005). 



Hnat et all (1200 3a). and the scaling in r is broken 



as we approach the characteristic substorm timescale of 1-2 hours. Remnants of the substorm 
signature will be present in the time series on timescales shorter than this. The behaviour of the 
peak of the PDF (P(y —►())) will also be dominated by uncertainties in the determination of the 
signal rather than its scaling properties. 

Figure 3 shows a plot of versus p for the AE time series in the same format as figure 2 
for the interval January 1978 to July 1979 comprising 7.5 x 10 5 samples. Plots of the structure 
functions used to construct figure 3 are shown in figure 4. The error bars on figure 3 are those of 
the best fit straight lines to Figure 4 rather than the possible range of straight line fits and as such 
are a minimal error estimate. 

We plot in figure 4(a) the raw result, that is d!6t and in figure 4(b) the conditioned approximation 
(flTt with Q = 20, the latter corresponding to the removal of less than 1 % of the data. From figure 
4 we see that no clear scaling emerges beyond the third order p = 3 until approximation dT7l is 
made. Clearly, if scaling is present, the ((p) obtained from the raw structure functions (figure 
4(a)) are not a good estimate. Once the data is conditioned, we find that Q = [10, 20] give almost 
identical estimates of ((p) which are weakly multifractal. For Q = 5 the are shifted slightly 
toward self similar scaling. The closeness of the conditioned results for the range Q = [5, 20], and 
their clear separation from the raw result, suggests that these are a reasonable approximate measure 
of the scaling properties of the time series. This procedure can be used to make quantitative 
comparisons between timeseries to this precision. Given the caveats above however, we cannot 
use this procedure to distinguish whether the time series is self affine or weakly multifractal, but 
can distinguish strong multifractality. 



4 Low Order Moments and Non Uniqueness: Comparison with a fractional Brownian sur- 
rogate. 

Equation <l 14l> expresses the PDF as an expansion in the moments to all orders. It follows that 
distinct timeseries can share the first few moments and therefore if scaling, may also share the 
same Hurst exponent and corresponding exponent of the power law power spectrum. Having 
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estimated the scaling exponent of the AE index as above we can construct a time series with the 
same second moment from a fractional Brownian motion to illustrate this. 

The fraction al Brownian walk was generated using the method described in Appendix 3 of 



Peters. 



( 1996). The algorithm takes a series of Gaussian random numbers and approximates a 
finite correlation time by weighting past values according to a power law function. In our case 
1024 Gaussian samples were used to create each increment of fractional walk. The resulting time 
series is comprised of 7.5 x 10 5 increments. 

Figure 5 shows the two time series, (i) the interval of AE analyzed above, and (ii) the fBm 
surrogate. The standard deviation versus r for the two time series is s hown in Figure 6. Th e 



power spectrum of AE (the raw, rather than the differenced variable)( c.f. usurutani et al. 



Takalo et al. 



(1990); 



1 19931) ). along with the er(r) and the structure functions, show a characteristic break 



on t imescales of 1-2 ho urs. O n times shorter tha n this, we can obtain a scaling collapse of the PDF 



(see 



HnatetaL 



l 2003ah . also 



HnatetaL 



(I20051) . Fluctuations on these timescales share the same 
second moment as the fBm. In Figure 7 we compare the PDF of these flu ctuations and we se e that 
these are very distinct; fBm is defined as having Gaussian increments AMandelbroti d2002h ') and 
this is revealed by the PDF whereas the AE increments are non-Gaussian. 

This is an illustration of the fact that the scaling in AE over this region is not necessarily due 
to time correlation, the "Joseph effect" for which Mandelbrot constructed fractional Brownian 
motion as a model. Indeed AE has almost uncorre cted differences at high frequencies, as indi- 



cated by its nearly Brownian / 1 power spectrum iTsurutani et al.\ dl99oL. Rather the scaling 



is synonymous with the heavy tailed PDF ("Noah effect") for which Mandelbrot (2002) earlier 
introduced a Levy model in economics. 

Finally, we plot in Figure 8 the ((p) versus p obtained from the structure function estimate dTTl 
with Q = 10 for both time series. We see from the plot that both time series are self affine and 
to within the uncertainty of the finite length time series, both share values of C{p) for the lowest 
orders in p. However the higher order structure functions reveal the distinct scaling of the two 
time series. 



5 Fokker-Planck model 

For completeness we now outline how the exponent a of a self affine time series leads to the 
functional form of P{y) via a Fokker- Planck model of the stochastic process x(t). Here we will 
consider an approach where scaling is achieved via transport coefficients that are functions of the 
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differenced variable y(t ). An alternative 



(y) coordinate ( see e.g. 



equivalent (e.g. 



Schertzer et aL 



approach is via fractiona 



(2001) 



Yannacopoulos and Rowlands 



Sh les ineer et al. 



( 1997k 



deriv atives for the dependent 



(1995|)). These are in principle 



We begin with a general form of the Fokker-Planck equation can be written 

dP 



dr 



V y (A(y)P + B(y)V y P), 



Gardinen dl986h : 
(18) 



where P = P(y, t) is a PDF for the differenced quantity y that varies with time r, A{y) is the 
friction coefficient and B{y) is related to a diffusion coefficient which we allow to vary with y. If 
we now impose the condition that solutions of (fT8l are invariant under the scaling given by (J9ji, 
then it is fo und that both A(y) a nd B(y) must have the form of power law dependence on y. Then 



as shown in 



Hnatetal 



(2003b), (HSJl takes the form: 



dP 



d_ 

dy 



y I y 



-l/a 



a P + b y 



dP 
dy 



(19) 



where ao and bo are constants, a is the scaling index derived from the data and P{y), y are 
unsealed PDF and fluctuations respectively, and where here we have explicitly insisted that the 
diffusion coefficient B{y) > 0. Importantly, in a physical system the scaling behaviour Q is 
expected to be strongly modified as y — ► 0, that is, at the peak of the PDF P(y) since for a 
sufficiently small difference between two measurements x(t), y(t,r) = x(t + r) — x(t) will be 
dominated by the uncertainties in those measurements. 

Written in this form equation dT9b immediately allows us to identity B(y) oc y 2 \ y \~ l l a 
and A(y) oc y \ y l^ 1 / . Solutions to d 19l > exist which are functions of y s = yr~ a only which 
correspond to stationary solutions with respect to r. We obtain these by the change of variables 
(P,y,r ^ P s ,y s ) of d: 

C | y s | « 



a 



bo dP s 

—y s -; — i- P a h — \y s \<* p s 

a dy s ao 



(20) 



This differential equation d2"0l can be solved analytically with a general solution of the form: 

C ' ~2 



«0 



X 



P s (y s ) 
y 3 | y' s I "o exp [ ^ | y 



exp 



b \y s \ a °/ b ° 

o?_ I / il/a 



a 
~~b~o 



\l/a 



d(y' s ) + k H(y s ), 



(21) 



'° I Vs \ a 

where ko is a constant and H(y s ) is the homogeneous solution: 

1 



H{5x s 



\ao/bo 



a i il/a 
~~b~o Vs 



(22) 
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Power law scaling for arbitrary y leads to singular behaviour of this solution at y — > 0. We do 
not however expect this to describe a physical system as y — > as discussed above. For JTil 
to describe a PDF we require that its integral is finite. We can discuss this by considering the 
behaviour close to the singularity: 



lim P(y s 



ao C 



y's I "° dy' s 
/2-i 



+ 



kn 



c + 



k 



a 
'■'11 



(23) 



Vs " I Vs \°° I Us 

The integral of d23t is finite for < ao/bo < 1 and < a < 1/2 (a subdiffusive process) so that 
within this range the integral of (12 1 b will be finite also as required. Outside of this range it can only 
be considered as an asymptotic solution. However, we can consider the generalization y — > y + e 
in the above, where e is a constant of magnitude that is small compared to, say, the values of <t(t) 
for the physical system under study. This eliminates the singular behaviour and corresponds (for 
y small) to the addition of low amplitude Gaussian noise as can be seen from the form of the 
corresponding Langevin equation d24l below. Physically this corresponds to a simple model for 
the statistical behaviour of the observational uncertainties in the data which may dominate as the 
differenced quantity y — > 0. 

Expression d2~TT i is then a family of solutions for the PDF of self affine time series. This provides 
a method to test for self affinity that does not directly rely on determining the scaling exponents 
to high order from the structure functions. Having determined the exponent a from the scaling 
of a low order moment (say, the standard deviation) one can then perform a scaling collapse on 
the PDF ; this should then a lso be described by the corresponding solution of (ITU (see 



aOOm - AHnat e.t al. 



Hnat et al. 



(2005) for examples). 



It is well know n that a Fokker Planck equation is simply related to a Langevin equation (see e.g. 



Gardinen &19&(b ). A nonlinear Langevin equation of the form 



d f t =P(y) + i(y)m, 



(24) 



where 0(y) is ay -d ependent force term and j(y) is a y -dependent noise strength, can be shown 



(Hna t et al. 



d2003bh ) to correspond to dT8l and in that sense to describe the time series. In d24l the 
random variable £(i) is assumed to be ^-correlated, i.e., 

<£(*)£(* + t) >=a 2 5(r). (25) 

Consistency with equation (0 is achieved in the data analysis by forming each time series y(t,r) 
with non-overlapping time intervals t. Defining Dq =< i 2 (t) > /2 we then obtain: 



iJo 



(26) 
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and 



a 



y\y\ 



(27) 



With a = 1/2 and clq = one recovers the Brownian random walk with (II 8b reduced to a diffusion 
equation with c onstant diffu sion coefficient. 

Interestingly, IBecIa (120011) has independently proposed a nonlinear Langevin equation where 
but not 7 varies with y. This yields leptokurtic PDFs of the Tsallis functional form. 

Finally the variable r in COD, and t in (Ell can be read in two ways: either as the renormalization 
variable of the stochastic variable y(t, r) or the time variable of x(t) since from $6^ t = 2 n A and 
with the choice N = 2 n we have x(t) = y n (t, r), r = t (n, N large). Thus d24t can be seen either 
as a prescription for generating a self- affine timeseries with scaling exponent a, or as describing 
the renormalization flow. 



6 Conclusions 

Empirical determination of the scaling properties and exponents of time series x(t) presents a 
formidable challenge in testing, and developing, a theoretical understanding of turbulence and 
other out-of-equilibrium phenomena. In this paper we have discussed the special case of self 
affme time series by treating the differenced variable y(t,r) = x(t + r) — x(t) as increments 
of a stochastic process (a generalized random walk). We have highlighted two complementary 
approaches to the data. 

The first of these is PDF rescaling; using a low order moment to determine a scaling exponent 
and then verifying whether this exponent collapses the PDFs of the differenced variable y(t, r) 
over the full range of y accessible from the data. As a corollary this collapsed PDF should also 
be well described by the solution of a Fokker-Planck equation which has power law transport 
coefficients. 

The second of these is using structure functions to determine the scaling properties of the higher 
order moments. In a finite length time series the higher order structure functions can be distorted 
by isolated, extreme events which are not well represented statistically. Using the example of a 
finite length Levy flight, we have demonstrated a method for conditioning the time series that can 
in principle recover the underlying self affine scaling. 

Finally, to highlight how both these methods are complementary in quantifying the scaling 
properties of the time series a fractional Brownian walk was constructed to share the same second 
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moment as an interval of the differenced AE index time series. The two timeseries were demon- 
strated to possess very different PDF of the differenced variable, and distinct structure functions. 
Both of these appro aches could in principle be generalized to multifractal time series (see e.g. 



Schertzer et al. 



<2001|)). 
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Figure Captions 



Fig. 1. Scaling collapse of the PDF of an fi = 1.8 Levy flight. 

Fig. 2. The effect of conditioning a Levy flight, -jc are the ((p) obtained from the raw time series, all other 
symbols refer to conditioned time series for different values of Q (see text). The conditioned results yield a 
scaling exponent a = 0.544 which corresponds to a Levy index of /i = 1.84. 

Fig. 3. Scaling exponents £(p) versus p for the AE index, shown in the same format as figure 2 

Fig. 4. Structure functions of the AE index estimated for orders p = [1, 6] by method (1161 (a) and by 
method O (b). 

Fig. 5. A ~ 1.5 year interval of AE data (upper trace) is shown alongside a surrogate fBm time series 
(lower trace) with the same second moment. The traces have been displaced for clarity. 

Fig. 6. Standard deviation of the time series differenced on timescale r plotted versus r for an interval of 
AE index data (see text) and an fBm time series constructed with the same second moment. The traces have 
been displaced for clarity. 

Fig. 7. PDF of the time series of AE, differenced on timescales less than one hour (□). The PDF of an fBm 
with the same second moment is shown for comparison (o). 

Fig. 8. Structure functions obtained by conditioning at 10 a for an interval of the AE index, and for a fBm 
constructed to share the same second moment. 
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Fig. 1. Scaling collapse of the PDF of an /i = 1.8 Levy flight. 
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Fig. 2. The effect of conditioning a Levy flight. + are the ((p) obtained from the raw time series, all other 
symbols refer to conditioned time series for different values of Q (see text). The conditioned results yield a 
scaling exponent a = 0.544 which corresponds to a Levy index of /i = 1.84. 
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Fig. 3. Scaling exponents £(p) versus p for the AE index, shown in the same format as figure 2 
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Fig. 5. A ~ 1.5 year interval of AE data (upper trace) is shown alongside a surrogate fBm time series 
(lower trace) with the same second moment. The traces have been displaced for clarity. 




Fig. 6. Standard deviation of the time series differenced on timescale r plotted versus r for an interval of 
AE index data (see text) and an fBm time series constructed with the same second moment. The traces have 
been displaced for clarity. 
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Fig. 7. PDF of the time series of AE, differenced on timescales less than one hour (□). The PDF of an ffim 
with the same second moment is shown for comparison (o). 
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Fig. 8. Structure functions obtained by conditioning at 10 a for an interval of the AE index, and for a fBm 
constructed to share the same second moment. 
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